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ABSTRACT 

We present high angular resolution (-0.94") ^^CO(l-O) Atacama Large MillimeterlSubmillimeter Array 
(ALMA) observations obtained during the 2014 long baseline campaign from the eruptive bipolar flow from the 
multiple XZ Tau stellar system discovered by the Hubble Space Telescope (HST). These observations reveal, 
for the first time, the kinematics of the molecular flow. The kinematics of the different ejections close to XZ Tau 
reveal a rotating and expanding structure with a southeast-northwest velocity gradient. The youngest eruptive 
bubbles unveiled in the optical HST images are inside of this molecular expanding structure. Additionally, 
we report a very compact and collimated bipolar outflow emanating from XZ Tau A, which indicates that the 
eruptive outflow is indeed originating from this object. The mass (3 x 10“^ M©) and energetics (Ekin = 3 x 
10^^ ergs) for the collimated outflow are comparable with those found in molecular outflows associated with 
young brown dwarfs. 

Subject headings: stars: pre-main sequence - ISM: jets and outflows - individual: (XZ Tauri) - individual: 

(Haro 6-15) 


1. INTRODUCTION 

XZ Tau is considered a classical T Tauri star (classified 
as a class I/II), which is composed of a binary system, XZ 
Tau (A/B). The binary system is separated by approximately 
0.3" (at a PA. of 154°), or 42 AU at a distance of 150 pc 
(Torres et al. 2009). The relative flux-ratio between the two 
components fluctuates at optical and radio wavelengths, prob¬ 
ably because of non-thermal (gyro)synchrotron emission as¬ 
sociated with magnetic activity (Krist et al. 2008). XZ Tau is 
exciting an optical outflow (HH 151) with a PA. = 15°, and 
which can be traced over about 10" on either side of the star 
as revealed by the studies of Mundt et al. (1990). 

Contrary to the long and well-collimated jet emanating 
from the nearby T Tauri star HL Tau (Krist et al. 2008; 
Lumbreras & Zapata 2014), XZ Tau is driving an outflow 
with a morphology reminiscent of numerous concentric ex¬ 
plosions or hot bubbles, as revealed by multiple observa¬ 
tions carried out with the Hubble Space Telescope (HST) 
(Krist etal. 1999, 2008). There are only two reported ad¬ 
ditional cases where this morphology {i.e. bubble-like) is 
present: the SVS 13 outflow (Hodapp & Chini 2014) and 
IRASI6293-2422(B) (Loinard et al. 2013). The nature of 
these bubble-like outflows (including that of XZ Tau) is still 
under debate. 

The expansion (tangential, or in the plane of the sky) veloc¬ 
ities of the elongated optical bubbles in the flow from XZ Tau 
are about 70-200 kms“^ (Krist et al. 2008), while the radial 
velocities (in the line-of-sight) are about 100-300 km s“^ ob¬ 
tained from optical images of forbidden lines of nitrogen [Nil 
46583] or sulfur [SII 46716] (Hirth et al. 1997). Comparing 
the radial velocities and proper motions of the outflow, it has 
been suggested that the bubbles are mostly in the plane of 
the sky with an angle of about 30° with respect to this plane 
(Krist et al. 2008). Multiple images from the HST have re¬ 
vealed that the eruptive outflow is episodic, probably excited 
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by velocity pulses in the outflow of XZ Tau A (Krist et al. 
2008). 

Krist et al. (2008) suggested that bubbles are ejected from 
XZ Tau A because this source is better aligned to the 
flow and XZ Tau B seems very faint (Hirth et al. 1997). 
Krist etal. (1999, 2008) proposed that possibly the perias- 
tron passage of XZ Tau B might produce the eruptive out¬ 
flow. However, the short dynamical ages for the bubbles 
cannot be explained by the very slow motions of the binary 
system XZ Tau (A/B) in the plane of the sky (Krist et al. 
2008; Carrasco-Gonzalez et al. 2009; Forganetal. 2014). 
Carrasco-Gonzalez et al. (2009) reported the detection of a 
close radio binary system in XZ Tau A with a separation of 
only 0.09" (13 AU), and probably in an elliptical orbit. They 
suggested that the periastron passage of one component, XZ 
Tau C (named S2 by Carrasco-Gonzalez et al. 2009), of the 
discovered close radio binary system could much better ex¬ 
plain the eruptions in the flow from XZ Tau. 

It is speculated that the mass of the binary system (A/C) 
is around I M©, with a period for the system of only 40 
yr (Carrasco-Gonzalez et al. 2009). Therefore, the A/C sys¬ 
tem was at its previous periastron in the 1980s, and will be 
again at this position around 2020. However, the data pre¬ 
sented in Carrasco-Gonzalez et al. (2009) are not sufficient 
to confirm or reject this hypothesis, so this explanation for 
the bubble-like outflow generation was tentative. In addition, 
Forgan et al. (2014) did not detect the XZ Tau C component 
previously reported by Carrasco-Gonzalez et al. (2009). More 
observations of XZ Tau A/C at around 2020 would be required 
to strongly test the periastron passage model for bubble-like 
outflow generation (Forgan et al. 2014). 

In this Letter, we present new, high angular resolution 
^^CO(I-O) line observations obtained with i\\Q Atacama Large 
MillimeterlSubmillimeter Array (ALMA) toward the eruptive 
outflow from XZ Tau. Our analysis reveals the kinematics of 
the eruptive flow, as well as a very compact bipolar jet ema¬ 
nating from XZ Tau A. This compact molecular outflow is at 
the base of the optical eruptive flow. 

2. OBSERVATIONS 
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Fig. 1.— Left: HST/ACS (narrow band FR656N, grey scale) image overlaid with the positions of the 3 mm continuum sources from the eruptive flow in 
the XZ Tau region. Right: 3 mm (contour and color image) ALMA images from the XZ Tau A/B system. The continuum ALMA images were taken from 
Partnership et al. (2015a). The red contours are 0.386, 0.772, 1.16 and 1.54 mJy beam“^ The synthesized beam of the ALMA continuum image is shown in the 
lower-left corner of the right image. The optical image has an offset of about 0.4" to the south in order to correct the proper motions since the epoch of the HST 
observations, and the World Coordinate System (WCS) uncertainties for the HST images. 


The CO observations were made with 32-35 antennas 
(Partnership et al. 2015a,b) of ALMA in 2014 October and 
November, during the ALMA science verification data pro¬ 
gram for the long baseline campaign. The purpose of this 
campaign was to verify and demonstrate ALMA’s perfor¬ 
mance on baseline lengths larger than 10 km. These observa¬ 
tions only included antennas with diameters of 12 meters. The 
independent baselines ranged in projected length from 12 to 
15,240 m. The phase center of these millimeter observations 
was at Q'y 2000.0 = 04‘'31'"38M26, dnoma = -18°13'57'.'04, the 
position of the nearby HL Tau. The primary beam of ALMA 
at 115 GHz has a FWHM of ~55 arcsec, so it included XZ 
Tau within the primary field of view. As XZ Tau was at the 
eastern side of our FWHM, we corrected the final continuum 
image for the primary beam attenuation. 

The ALMA digital correlator was configured to observe the 
^^CO(l-O) line at a rest frequency of 115.27120 GHz. The 
spectral window to observe this line was configured with a 
total bandwidth of 234375.0 kHz and was divided into 3840 
channels, providing a spectral channel resolution of 61.035 
kHz or 0.15 km s“^ However, the spectral resolution in the 
data cube was smoothed to 0.25 km s“\ in order to have a 
better sensitivity in each velocity channel. 

Observations of the quasar J0510-f 1800 provided the abso¬ 
lute scale for the fiux density calibration, while the quasar 
J0431-f2037 provided the gain calibration. The quasar 
J0510-f 1800 was also used for the bandpass calibration. 

The data were calibrated, imaged, and analyzed using 
the Common Astronomy Software Applications (CASA), see 
Partnership et al. (2015a,b). To improve the surface bright¬ 
ness sensitivity, the ^^CO(l-O) data were tapered to a reso¬ 
lution of about 1.0". The resulting r.m.s. noise for the line 
emission was about 10 mJy beam“^ per velocity channel at an 


angular resolution 0.98" x 0.90" with a PA. of -8.66°. We 
used a robust parameter of 0.5 in the CLEAN task, in order to 
obtain a better signal-to-noise ratio, at the cost of losing some 
spatial resolution. 

The HST/ACS (narrow band FR656N) image from XZ Tau 
was obtained from the HST archive and this is only used as 
a reference for the ALMA data. This HST image has been 
extensively discussed in Krist et al. (2008). 

3. RESULTS AND DISCUSSION 

In Figure 1 (right panel), we present the resulting 3 mm 
continuum ALMA image of the XZ Tau. The positions of 
both the XZ Tau A and XZ Tau B components are shown 
in the left panel of Figure 1, superposed on an optical 
HST/ACS (narrow band FR656N) image of the eruptive out- 
fiow (Krist et al. 2008). This figure reveals the millimeter 
thermal emission arising from the binary system (A/B) lo¬ 
cated in XZ Tau and locates the continuum sources with re¬ 
spect to the outfiow. The continuum emission at 3.3 mm is 
presented and discussed in Partnership et al. (2015a). 

We present the first moment (that is, the velocity-weighted 
intensity) map of the CO(l-O) line emission in Figure 2 (the 
velocity range of integration is from -2.25 to -r 11.75 km s“^). 
This image reveals the gas kinematics of the eruptive outfiow 
associated with XZ Tau. We have additionally overlaid the 
CO velocity field on the HST/ACS optical image (using the 
spectral filter FR656N) and we have indicated with points the 
positions of the ALMA 3.3 mm continuum sources. This im¬ 
age shows the structure of the CO gas in the entire outfiow and 
in the vicinity of the optical jet. These observations there¬ 
fore reveal, for the first time, the kinematics of the molecu¬ 
lar flow. The systemic LSR velocity of the cloud where XZ 
Tau is located is ~6.8 km s“^ (as derived from the ^^CO ob- 
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Fig. 2.— ALMA first moment map (or the intensity weighted velocity map) 
of the CO(l-O) line emission from the bubble-like outflow in XZ Tau overlaid 
in grey contours with an HST/ACS optical image (using the spectral filter 
FR656N). We have marked with red dots the positions of the ALMA sources. 
The grey contours range from 0.05% to 0.2% in steps of 0.06% of the peak of 
the emission. The synthesized beam of the ALMA CO(l-O) image is shown 
in the lower-left corner. The optical image has an offset of about 0.4" to the 
south in order to correct for the proper motions since the epoch of the HST 
observations. The dashed arrows mark the position of the cuts (C1/C2) and 
their orientations toward the positive values (see Fig. 3). 


servations of Calvet et al, 1983, see also Partnership et al. 
2015a). The velocity centroid of the expanding structure is 
~5.8 km s“^ (see Fig. 3). Then, as a whole, this structure 
is slightly blueshifted (by about 1 km s“^) with respect to 
the ambient gas. There is a clear velocity gradient with a 
southeast-northwest orientation close to the position of XZ 
Tau. This molecular gas could be associated with expanding 
gas observed in the hot optical bubbles revealed by the HST 
at optical wavelengths. The CO gas associated with the older 
eruptions (located far from the XZ Tau multiple system) of 
the flow is mostly associated with low velocity gas close to 
the ambient material, presumably because it is decelerating. 

Figure 3 shows the position-velocity (PV) diagrams for two 
cuts (Cl and C2) across (PA. = 124°) and along (PA. = 
24°) the eruptive outflow (Figures 2 and 3), which reveal 
the kinematics of the molecular gas with respect to the po¬ 
sition. The intersection point of these position-velocity dia¬ 
grams is placed at the estimated centroid of the low-velocity 
emission, about 1 arcsec to the north of the stellar system. 
These position-velocity cuts reveal an expanding and rotating 
molecular structure. The expansion is evident in the velocity 
components at about -\-5 and -f7 km s“^ and the rotation in 
the (south)east-(north)west velocity gradient. Deeper obser¬ 
vations are needed to reflne our kinematic interpretation for 
this rotating and expanding structure. 

Three more signiflcant features in the PV cuts are present: 
a central feature at V ^ 6 km s“^ and 0 positional offset with 
respect to the intersection of the cuts, extended emission close 
to (most notable in the C2 cut), and fainter higher veloc¬ 
ity blueshifted and redshifted emission mainly along the C2 
cut. This higher velocity is discussed below. 

Finally, in Figure 4 additionally to the 3 mm continuum and 
the optical HST images, we overlaid the integrated intensity 
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Fig. 3.— Position-velocity diagrams across (P.A. = 124°) and along (PA. 
= 24°) the molecular outflow in XZ Tau: radial velocity as function of on- 
the-sky distance. In Figure 2, we show the position of both cuts (Cl and C2). 
The angular offsets are in arcseconds and the LSR radial velocities in km s“^. 
In both panels the contours range from 40% to 90% of the peak emission, 
in steps of 10%. The peak of the line emission is 0.11 Jy beam“k The 
synthesized beam is 0.98" x 0.90" with a P.A. of -8.66°, and the spectral 
resolution is 0.25 km s“k 


(moment 0) to the higher velocity map of the ^^CO(l-O) emis¬ 
sion. In the image, the blue and red colors correspond to the 
moderately-high velocity blueshifted and redshifted gas emis¬ 
sion, respectively. We integrated LSR velocities from -fO.75 
to -f4.5 km s“^ for the blueshifted lobe, and from -h8.25 to 
-f 10.75 km s“^ for the redshifted one. This image reveals a 
compact bipolar molecular outflow emanating from XZ Tau. 
This outflow is narrow and collimated, similar in morphol¬ 
ogy to a high-velocity molecular outflow (Zapata et al. 2005). 
However, the molecular gas seen at these moderate velocities 
is likely not the original gas ejected from the star/disk system 
but mostly entrained ambient material. The jet shows a wavy 
morphology, observed mostly in the blueshifted lobe (see Fig¬ 
ure 4), possibly caused by a binary associated with component 
A (Carrasco-Gonzalez et al. 2009). The P.A. of the jet is ap¬ 
proximately 20°. 

Morphologically, the blueshifted lobe seems to connect 
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Fig. 4. — ALMA high-velocity integrated CO(l-O) emission (zero moment map) overlaid on an HST/ACS optical image (narrow band FR656N) from the 
XZ Tan region. The white contours mark the 3 mm continuum emission from the binary in XZ Tau. The blue contours represent molecular emission that is 
approaching to us, while the red contours represent the emission that is receding. The integrated range of velocities for the blue contours is from -I-0.75 to -I-4.5 
km s“^ and for the red emission is from -I-8.25 to -i-10.75 km s“8 The white contours are in percent of the peak emission starting from 10% to 90% in steps of 
10%. The intensity peak emission is 1.8 mJy beam“8 The blue and red contours are similar, starting from 30% to 90% in steps of 10% for the blue contours 
and 15% to 45% in steps of 5% for the red ones. The peaks of the emission are 0.052 (blue) and 0.028 (red) Jy beam“^ km s“8 The synthesized beams of the 
ALMA CO(l-O) (blue contour) and 3 mm continuum (black contour) images are shown in the lower-left corner. The optical image has an offset of about 0.4" to 
the south in order to correct for the proper motions since the HST epoch, and the WCS uncertainties for the HST images. 


with one optical knot far from the origin. The redshifted side 
is more compact. The molecular outflow is ejected from the 
XZ Tau A object, as first suggested by Krist et al. (2008) us¬ 
ing HST optical observations. These ALMA data confirm the 
suggestion of Krist et al. that the origin of the flow is XZ 
Tau A. It is interesting to mention that the orientations of the 
CO outflow here reported and the optical jet/outflow as seen 
in [SII] and [Nil] lines (Mundt et al. 1990, Hirth et al. 1997) 
agree, but the blue and redshifted sides appear to be inverted. 
This could suggest that they are independent features in the 
XZ Tau system. Another possible explanation for this discrep¬ 
ancy is that while the optical jet traces fast, collimated, and 


redshifted atomic gas, the CO is tracing part of an entrained 
molecular flow that, if wide enough, could have a blueshifted 
component (Takakuwa & Kamazaki 2011). 

Assuming that the ^^CO(l-O) line emission is optically 
thin and in Local Thermodynamic Equilibrium, we esti¬ 
mate the outflow mass using the following equation (e.g., 
Scoville et al. 1986; Palau et al. 2007), which for the (1-0) 
transition is given by: 


Mh2 


= 3.04xl0“^^rex^^Xi/2 
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where we used 2.8 as mean molecular weight, and is the 

CO 

abundance ratio between the molecular hydrogen and the car¬ 
bon monoxide (-10^^, e.g., Scoville et al. 1986), Tex is in units 
of K and assumed to be 10 K, J lydv is the average intensity 
integrated over velocity, ^maj and ^min are the projected major 
and minor axes of the outflow lobe, and D is the distance to the 
source (150 pc). We then estimate a mass for the molecular 
collimated outflow powered by XZ Tau A of approximately 3 
X 10“^ Mo. 

This value is quite low compared to the molecular mass 
estimated for outflows excited by young low-mass proto¬ 
stars and observed with interferometers: HH 211, Palau et al. 
(2006); HL Tau Lumbreras & Zapata (2014); B5-IRS1 
(Zapata et al. 2014); and DG Tau B (Zapata et al. 2015), and 
even lower, but comparable, to those observed in molecular 
outflows from sub-stellar objects (e.g. Phan-Bao et al. 2014). 
We also estimated a kinetic energy of 3 x 10^^ ergs. 

Taking a radial velocity for the CO of 5 km s“^ (which are 
the velocities for the moderately-high velocity collimated out¬ 
flow), a size of 1" (see Figure 4) and an angle with respect to 
the plane of the sky of 30°, we estimate a dynamical age of 
about tdyn = (size of the outflow in the plane of the sky / radial 
velocity) x tan[^] = (2.3x10^^ cm / 5 xlO^ cm s“^) x tan[^] 
= 148 tan[^] years for the collimated molecular outflow. In 
this equation 6 is the angle of the outflow with respect to the 
plane of the sky. Using the determinations and respective er¬ 
rors for the radial and transversal velocities (Hirth et al. 1997; 
Krist et al. 2008) for the optical jet, we And that 6 can be in 
a wide range of values, from 13 to 61 degrees. This implies 
an age estimate between 33 and 266 years, with the smaller 
value close to the estimated age of 34 years for the bubble by 
the time the ALMA data was taken. Then, the ejection of the 
optical bubbles and the creation of the molecular collimated 
outflow could have been simultaneous, however the fact that 
the orientations are inverted (for the collimated outflow and 
the bubbles) suggests they are probably different structures. 

4. CONCLUSIONS 

We have analyzed sensitive CO line ALMA observations in 
the millimeter regime from the young multiple star XZ Tau. 
Our conclusions are as follows: 

• The ^^CO(l-O) line observations show for the first time 
the kinematics of the molecular outflow. The kinemat¬ 
ics of the different ejections close to XZ Tau reveal a 
slow (AV = 2 km s“^) and expanding/rotating structure. 
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This structure appears to be related with the multiple 
eruptive bubbles unveiled in the optical HST images. 

• The observations also reveal the presence of a compact 
and collimated bipolar molecular outflow (reminiscent 
of a jet) with a northeast-southwest orientation (P.A. ~ 
20°) emanating from XZ Tau A. This outflow is prob¬ 
ably tracing the base of the optical eruptive flow, with 
a similar orientation. The mass (3 x 10“^ M©) and en¬ 
ergetics (3 X 10^^ ergs) for the molecular collimated 
outflow are smaller than those values found in other 
molecular collimated outflows associated with young 
low-mass stars. 

• The observations confirm the origin of the eruptive out¬ 
flow as XZ Tau A. These observations argue in fa¬ 
vor of the model presented by Carrasco-Gonzalez et al. 
(2009) with regard to the periodic nature of the outflow, 
however, as discussed in Forgan et al. (2014), addi¬ 
tional observations at different wavelengths are needed, 
especially dedicated observations to follow up the mo¬ 
tions, and the outflow ejections from the components in 
the XZ Tau multiple star system. This will help to iden¬ 
tify firmly the physical mechanisms behind the ejection 
of the bubble-like flow. 
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